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\Q • Abstract. We report the first observation of the reactions Au + Au — > Au + Au + p° 

and Au+Au^ Au* +Au* + p° with the STAR detector. The p are produced at small 
• perpendicular momentum, as expected if they couple coherently to both nuclei. We 

discuss models of vector meson production and the correlation with nuclear breakup, 
and present a fundamental test of quantum mechanics that is possible with the system. 
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1. Introduction 

With their large charge (Z), relativistic heavy ions can interact electromagnetically, even at 
impact parameters b much larger than twice the nuclear radius In these ultra-peripheral 
collisions, the two nuclei act as sources of fields. The electromagnetic fields have a long 
range, so two-photon and photon-Pomeron interactions have large cross sections. Exclu- 
sive Pomeron-Pomeron reactions have a small cross section because of the short Pomeron 
range. The amplitudes for the photon or Pomeron emission from each nucleon add. The 
momentum transfer is small enough that the amplitudes all have the same phase, so the 
coupling is coherent over the entire nucleus. The cross section for two-photon interac- 
tions goes as Z 4 (4 x 10 7 with gold ions). For photon-Pomeron interactions, o ~ Z 2 A 2 for 
'heavy' states like J /\\f and o ~ Z 2 A 4 / 3 for lighter mesons. This scaling leads to large cross 
sections. The coherence constrains the fields to have perpendicular momentum pj <H/Ra 
and longitudinal momentum pn < yfi/R^. Final states can have masses up to 2y/z/^ (6 
GeV at RHIC). The final state p T is less than — lhjR A . 

Photon-Pomeron or photon-meson interactions produce vector mesons, as in Fig. 1(a). 
The Pomeron is a colorless object with the same quantum numbers as the vacuum. Photon- 
Pomeron interactions can be described in terms of elastic scattering. A photon from the 
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electromagnetic field of one nucleus fluctuates to a qq state. This state then elastically scat- 
ters from the other nucleus, emerging as a vector meson. Vector meson dominance allows 
us to calculate the rate for this process by treating light qq states directly as vector mesons. 
At RHIC, p°, CO, (j) and are copiously produced, and the \|/' should be observable. The 
upsilon family may be observable with lighter ion beams. 



Fig. 1. The Feynman diagrams for (a) p° production, and (b) p° production with nuclear 
excitation. Although other diagrams can lead to p° production with nuclear excitation, (b) 
is believed to be dominant. The dashed line shows the factorization following Eq. (2). The 
number of photon lines in (b) is not fixed. Two photons are shown to demonstrate how the 
excitations of the two nuclei are independent. 

Two-photon interactions include e + e~ pair production, single meson production, and 
meson pair production. Since the coupling Za (0.6 for gold) is large, e + e~ pair production 
is an important probe of quantum electrodynamics in strong fields. Coupling to 2-photons 
is sensitive to a mesons internal charge (quark or gluon content). These physics topics are 
reviewed elsewhere[ 1 ] [2] . 

These reactions can be studied experimentally by selecting events with low multiplic- 
ities and small total p T . We present a study of p° production with and without nuclear 
excitation, based on data taken with the STAR (Solenoidal Tracker at RHIC) detector at 
RHIC at a center of mass energy of \/Snn = 130 GeV per nucleon pair. STAR is a collab- 
oration of about 320 physicists from 31 institutions in 8 countries. 

2. Expected Rates 

The cross section to produce a meson V is the convolution of the photon spectrum from one 
nucleus witho(yA — ► VA)[3]. The photon flux from one nucleus is given by the Weizsacker- 
Williams virtual photon approach; the small photon virtuality is negligible here. For a 
photon-Pomeron or two-photon interaction to be visible, there cannot be an accompanying 
hadronic interaction. This requirement is similar (but not identical) to requiring b > 2Ra; 
it reduces the photon flux by roughly a factor of 2. 




a) 



b) 
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The photon-nucleus cross section can be determined by a Glauber calculation that 
uses the a(yp — > Vp) as input. The rate depends on how vector mesons interact with 
nucleons. The RHIC vector meson production rates can be used to determine the vector 
meson-nucleon cross section. The cross sections are very large, 380 mb for p° production 
with gold ions at V^aw = 130 GeV. This is about 5% of the hadronic cross section. At 
design energy and luminosity, RHIC is expected to produce 120 p° each second. The 
production rates are high enough that rare meson decays should be accessible, and excited 
vector mesons such as the p*, co* and (j)* can be studied. 

Photons can also fluctuate to virtual 7C + 7C~ pairs. One of the Jt can interact with the 
target nucleus, and the pair can become real. The amplitude for y — > k + k~ is independent 
of the pair invariant mass m XK . Direct 7t + 7t~ production interferes with p° production, so 
their amplitudes add[4]; the cross section is 
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where A and B are the amplitudes for p and TOT production. The p width T p must be 
corrected for the increasing phase space as m %% increases. Several methods have been 
proposed; we take F p = Fq(p* / Po) 3 M p /M nn , where p* is the pion momentum in the jot 
rest frame, and p* Q = p* when M nn = M p [5] and M p = 768 MeV/c 2 , T p = 151 MeV/c 2 [6]. 
The 7C+7C~ spectrum also includes a small component from co — > 7t + 7t~ which is neglected 
here. 

As Eq. (1) shows, the p and TOT interference is constructive for m nn < m p . Around 
m KK = mp, the p component shifts phase by 180°, and at higher masses, the interference is 
destructive. This skews the overall m KK distribution. 



2.1. Vector Meson production with Nuclear Excitation 

Vector meson production can be accompanied by nuclear excitation, as in Fig. 1(b). For 
this, diagrams like Fig. 1(b) are expected to be dominant. One or more additional ex- 
changed photons can excite one or both nuclei to a giant dipole resonance. Higher collec- 
tive excitations or nuclear breakup are also possible. However, photons directly involved 
in two-photon or photon-Pomeron interactions are unlikely to also excite the emitting nu- 
cleus [7]. Pomeron emission (elastic scattering) is less likely to cause collective excitation 
than photon emission because Pomerons couple identically to protons and neutrons. Ne- 
glecting other production diagrams, the production probability factorizes and 

= J d 2 bP p (b)P 2GD R(b)[l-P H AD(b)] (2) 

where P p {b) is the b dependent probability of production a p, around 0.5% for b = 2Ra- 
The b dependence comes entirely from the photon flux. 

The single nucleus excitation probability is PcDR(b) = 1 — exp (— S/b 2 )[l], where S 
includes terms for different types of photo-disintegration. GDR excitation is the largest, but 
higher excitations and other photonuclear reactions also contribute. For gold at RHIC, S ~ 
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150 fm 2 . The probability of exciting both nuclei, P2GDit(b) ~ PcDR(b 2 ), with Pigdr^ = 
2Ra) ~ 30%[8]. These relationships assume that the excitation of the two nuclei occurs 
independently. 

Phad is the hadronic interaction probability, which can be approximated as 1 for b < 
2Ra\ otherwise. The cross section for vector meson production with nuclear excitation is 
expected to be more than an order of magnitude smaller than for exclusive vector meson 
production. 

2.2. Interference and the Vector Meson pr Spectrum 

Photons have a long range, and Pomerons (elastic scattering) have a short range, so vec- 
tor meson production must occur inside or near (within 1 fermi) one of the nuclei. Either 
nucleus can serve as the emitter or scatterer; the two possibilities are indistinguishable, so 
their amplitudes add and the two sources interfere. At mid-rapidity (y = 0), the contribu- 
tions from the two sources are equal. 

Vector mesons have negative parity, so the cross section at y = is [9] 

c(p T ,y = 0,b) = 2A 2 ( PT ,y = 0,b)(l - cos[p-b}). (3) 

where A is the production amplitude, b the impact parameter vector, and p the vector meson 
momentum. For a given b, a oscillates with period Apr =h/b. For p- b <C 1, the interfer- 
ence is destructive, and o is small. Of course, b is not observable in heavy ion collisions, 
so Eq. (3) must be integrated over all b. When this is done, the oscillations wash out ex- 
cept for pt < (b), where the cross section drops sharply. Fig. 2 shows the expected cross 
section as a function of pj. 



CM 




p ± [GeV/c] 

Fig. 2. Expected p° pj spectrum for gold on gold collisions at 200 GeV per nucleon, with 
(solid line) and without (dotted line) interference. 

This interference is of special interest because vector mesons decay quickly, before 
travelling the distance b required so their wave functions can overlap. The decay distance 
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yPcx is less than 1 fermi, far less than the (b) <~ 40 fermi at RHIC (and (b) <~ 300 fermi 
at the LHC). So, the vector mesons cannot interfere. However, their decay products can 
interfere, as long as the system wave function retains information on all possible decay 
amplitudes long after the decay occurs. Otherwise, for example, a J/\\t — > e + e~ decay 
at one nucleus couldn't interfere with a J/\y — ► hadrons at the other. The observation of 
interference will be a clear demonstration that the wave function doesn't collapse until 
after it is observed. Here, the 'observation' occurs via an interaction with the beampipe 
surrounding the collision region. 

3. Data Collection, Triggering and Analysis 

We have studied exclusive p° — > 7t + 7t~ production with the STAR detector. Data on gold- 
on-gold collisions at \/Snn = 130 GeV was collected during the Summer 2000 run. In this 
data, collisions could occur every 210 nsec. At STAR, the luminosity reached a maximum 
of about 2 x 10 25 /cm 2 /s. 

STAR detects charged particles in a 4. 2 meter long time proj ection chamber (TPC) [10]. 
The TPC has an inner radius of 50 cm, and an outer radius of 2 m. The TPC is centered 
around the interaction region. The pseudorapidity acceptance for charged particles depends 
on the particle production point; for the data discussed here, the interactions were spread 
longitudinally (in z), with a z = 90 cm. The pseudorapidity (y) acceptance of the trigger 
and track reconstruction depend on the individual vertex position. For a vertex at z = 0, 
the trigger was sensitive to charged particles with |y| < 1; the tracking covered a somewhat 
larger solid angle. A solenoidal magnet surrounds the TPC; for this data, the field was 
0.25 T. In this field, the TPC momentum resolution was about Ap/p = 2%. Tracks with 
Pt > 100 MeV/c were reconstructed with good efficiency. Tracks were identified by their 
energy loss in the TPC; this dE/dx was measured with a resolution of about 8%. 

The TPC is surrounded by a cylindrical central trigger barrel (CTB). For tracks from 
the center of the interaction region, it is sensitive to tracks with |y| < 1.0. This barrel 
consists of 240 scintillator slats, each covering Ay = 0.5 by A(|) = Jt/30. The scintillator 
light output is digitized to 8 bits of accuracy on each crossing. In the 0.25 T magnetic field, 
the scintillator was sensitive to charged particles with pr > 130 MeV/c. 

Two zero degree calorimeters (ZDCs) at z = ± 1 8 meters from the interaction point de- 
tect neutrons from nuclear breakup [11]. These calorimeters are sensitive to single neutrons, 
and have an acceptance of close to 100% for neutrons from nuclear breakup. 

The trigger hardware has several levels. The initial Level decision uses lookup tables 
and field programmable gate arrays to initiate TPC readout about 1 .5 ps after the collision. 
The other level used here, Level 3, is based on on-line reconstruction using a small farm of 
processors [12]; a level 3 acceptance triggered event building and data transfer to tape. 

3.1. Triggering and Data Collection 

We studied p° production with two separate triggers. The topology trigger was designed to 
trigger on p° decay products detected in the CTB system. A minimum bias trigger used the 
ZDCs to select events where both nuclei dissociated. 
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The topology trigger was sensitive to a number of different backgrounds. The major 
ones were cosmic rays, beam gas interactions, and debris from upstream interactions [13]. 
The latter refers to beam gas and/or beam-beampipe collisions far upstream from STAR; 
the usual manifestation was one or more tracks in the TPC roughly parallel to the beampipe, 
often accompanied by softer debris, perhaps from interactions in the TPC wheel or mag- 
netic pole tip. 



/ Top Veto 




\ Bottom Veto 
Cosmic Ray 



Fig. 3. Schematic view of the topology trigger. 

Fig. 3 shows how the topology trigger divided the CTB into 4 azimuthal quadrants. It 
selected events with at least one hit in the North and South sectors. The top and bottom 
quadrants were used as vetos to reject most cosmic rays. The trigger rate varied from 
20 to 40 Hz, depending on luminosity. Triggered events were scrutinized by the Level 
3 trigger, which reconstructed charged tracks online. Events with more than 15 tracks, 
or with a vertex far outside the interaction diamond were rejected. The multiplicity cut 
removed central collision events which left remnant energy in the trigger detectors; this 
remnant energy could cause a topology trigger on subsequent beam crossings. The vertex 
cut removed debris from upstream interactions, cosmic rays, and beam gas events. These 
cuts rejected 95% of the events, leaving 1-2 events/second for event building and storage. 
In about 7 hours of data taking, we collected about 30,000 events with this trigger. The 
luminosity ranged from about 2 x 10 24 /cm 2 /s up to (briefly) 2 x 10 25 /cm 2 /s, 

The minimum bias trigger required a coincidence between the two zero degree calorime- 
ters. The thresholds were set so that the efficiency was high for single (or more) neutron 
deposition. This trigger ran throughout the STAR data taking. This analysis is based on 
about 400,000 events. 

3.2. Exclusive p° Analysis 

Our analysis selected events with exactly two reconstructed tracks in the TPC. The tracks 
were vertexed with a low-multiplicity vertex finder, and required to form a vertex within 
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2 cm of the TPC center in x and y, and within 2 m in z. The vertexer projected the tracks 
back from the TPC to the interaction region with allowance for multiple scattering. Events 
were accepted if the two tracks were consistent with coming from a single vertex. 

As the \z\ position of the vertex increases, the solid angle covered by the TPC de- 
creases, so events at larger |z| are subject to more background by higher-multiplicity pro- 
cesses where one or more tracks is missed. The major analysis backgrounds were grazing 
nuclear collisions, incoherent photonuclear events, and beam gas events. Some cosmic rays 
also remained. 

To reject the remaining cosmic rays, track pairs were required to be at least slightly 
acoplanar, with a 3-dimensional opening angle less than 3 radians. This necessary cut 
reduced the p° reconstruction efficiency near y = 0, where the two pion tracks are nearly 
back-to-back. 

Track dE j dx was required to be consistent with that expected for a pion. Unfortu- 
nately, in the kinematic range for p — * TOT, the K and e dE / dx bands overlap, and only a 
few events were rejected. 

Fig. 4(a) shows the pj spectrum of the topology triggered pairs that pass these cuts, 
for unlike sign (dots - net charge 0) and like sign pairs (histograms). For the unlike pairs, a 
large peak is visible at pj < 100 MeV/c. This is consistent with production that is coherent 
with both nuclei; the events with pj < 100 MeV/c are considered our signal. The like sign 
pairs have no such enhancement, and can serve as a background sample. The like-sign pairs 
have been normalized to match the unlike sign in the signal-free region 0. 1 GeV < m Kll < 
1.0 GeV; this entailed scaling them up by a factor of 2.1. 



Transverse Momenlum | Invariant Mass lor pt<0.1 




pT (GeV) M(n+it-) (GeV) 



Fig. 4. (a) The pj spectrum of topology triggered 2-track events, (b) The m%% spectrum of 
2-track events with pj < 100 MeV/c. The points are oppositely charged pairs, while the 
histograms are the like-sign background, scaled up by a factor of 2.1. 



Fig. 4(b) shows the invariant mass of the pairs with pj < 100 MeV/c. The points 
are the unlike sign events, while the hatched histogram are the scaled like-sign pairs. The 
like-sign pairs are concentrated at relatively low masses, while the net charge pairs have 
a peak around the rho mass. We will consider the peak shape later, for the summed m nn 
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spectrum from both analyses. 

3.3. p° with Nuclear Excitation Analysis 

Events were selected by the minimum bias trigger, which took data throughout the sum- 
mer. The acceptance for these events was independent of the p° kinematics and the CTB 
calibration, simplifying the interpretation. 

The event selection was the same as for the topology trigger sample. Fig. 5(a) shows 
the pr distribution of 2-track events. The opposite sign pairs show the same peak at pj < 
100 MeV/c as the topology triggered data. The like sign background shows no peak, and, 
in fact, goes to zero as pj — > 0, as expected for a phase-space distribution. The like-sign 
pairs were scaled up to match the un-like sign pairs in the region 0.1 GeV < < 1.0 
GeV; this scaling factor was 2.3. 



Transverse Momentum | Invariant Mass for pT<0.1 




pT(GeV) M(jl+jr-) (GeV) 



Fig. 5. (a) The pj spectrum of minimum bias 2-track events, (b) The m KK spectrum of 
2-track events with pj < 100 MeV/c. The points are oppositely charged pairs, while the 
histograms are the like-sign background, scaled up by 2.3. 

This data contains little background from beam gas and incoherent photonuclear events, 
which are unlikely to deposit energy in both ZDCs. However, the background from grazing 
nuclear collisions is larger. 

Fig. 5(b) shows the m Kn distribution for pj < 100 MeV/c in the same data, along with 
the scaled like-sign background. The signal is peaked around the p mass. 

To improve the statistics, the minimum-bias and topology triggered samples are com- 
bined into a single m KK spectrum, shown in Fig. 6. The spectrum is fit to Eq. (1). In 
the fit, Mp and T p are fixed to minimize the number of free parameters. However, if the 
mass and width are allowed to float, the found values are consistent with the actual values. 
We find the relative size of the %% contribution, \B/A\ = 0.89 ± 0.07/VGeV (statistical 
error only). This ratio is somewhat larger than that found by the ZEUS collaboration, 
\B/A \ — 0.81 ±0.04/\/GeV at the same momentum transferf 14]. This is unexpected, since 
the TOT component should have a higher nuclear absorption cross section, so \B/A\ should 
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Fig. 6. The combined mjnt invariant mass spectrum, with a fit to p and nit components. The 
histogram is the data. The Breit-Wigner dashed curve is the p, the flat dashed line direct 
7t + 7t~ , and the solid curve the combined fit (including the interference). 



drop as the target atomic number increases [15]. 

The background (not yet subtracted) may include a component from electromagnetic 
e + e~ pair production, which would also be produced with low pj. The background is 
concentrated at low m nn , so a subtraction will reduce the apparent TOT component. 

Fig. 7 compares the ZDC values from the two datasets. The minimum bias data shows 
a clear single-neutron peak, with a significant component at higher energies. In contrast, 
the topology triggered data shows almost no energy deposition. This shows that the two 
processes are clearly distinguishable. 

4. Conclusions 

We have observed for the first time the reactions Au+Au^>Au+Au + p° and Au+Au^> 
Au*+Au* + p°. The p are produced at small perpendicular momentum, showing their 
coherent coupling to both nuclei. 

In the coming year, we expect to greatly expand the physics reach of the STAR ultra- 
peripheral collisions program. The improvement will come from many factors, among 
them the increasing capability of the STAR trigger, the additional detectors being installed 
in STAR, the increased beam energy, and the increased luminosity and running time. We 
plan to run the peripheral collisions trigger in parallel with the central collisions trigger(s), 
greatly increasing our data collection capabilities. The increased statistics will allow us to 
definitively detector or rule out interference among the two vector meson production sites. 
We will also be able to study two-photon production of mesons, and photoproduction of 
excited vector mesons. 
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Fig. 7. (a) shows the ZDC energy distribution (arbitrary units) for the topology triggered 
data. The peak at 4 is the pedestal, (b) shows the ZDC energy distribution for minimum bias 
data. The peak at 10 corresponds to single neutrons, while the higher-energy component 
corresponds to multiple neutron emission. 



Notes 



a. The collaboration membership is listed at http://www.star.bnl.gov/STAR/smdJ/collab/collab 
_sci_apr01.ps. 
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